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We investigate strongly correlated spin-1 ultracold bosons with antiferromagnetic interactions 
in a cubic optical lattice, based on bosonic dynamical mean-field theory. Rich phase diagrams of 
the system are mapped out at both zero and finite temperature, and in particular the existence 
of a spin-singlet condensate is established. Interestingly, at finite temperature, we find that the 
superfluid can be heated into a Mott insulator with even (odd) filling via a first- (second-) order 
phase transition, analogous to the Pomeranchuk effect in 3 He. Moreover, for typical experimental 
setups, we estimate the critical temperatures for different ordered phases and our results suggest 
that direct experimental observation of these phases is promising. 
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Spin correlations play an important role in understand¬ 
ing ground states of spinor bosons [U, and may lead 
to collective phenomena that are absent in the scalar 
condensate [2|. In recent years, spinor ultracold bosons 
have been successfully realized in experiments Q, and 
widely studied including spin mixing [i| , spin waves Q, 
spin dynamics Q, spin textures r] and phase transi¬ 
tions |8j. Experimental realizations open a new path 
to studying quantum phenomena of spin-correlated con¬ 
densates and competing insulating states, and also mo¬ 
tivate theoretical studies on spinor Bose gases in an op¬ 
tical lattice [S£l3. While quantum Monte-Carlo simu¬ 
lations and DMRG were available in one spatial dimen¬ 
sion [gi], up to now, these strongly correlated states in a 
three-dimensional (3D) optical lattice have been mainly 
investigated via lattice-gauge-theory predictions 101, ef¬ 
fective spin models in the strong-coupling limit 12, [h| 
or static mean-field theory |l4Hlq. Beyond static mean- 
field theory, a complete phase diagram of spinor bosons 
in a 3D optical lattice for the full range from small to 
large couplings is still unknown. Further open issues 
are the existence of a spin-singlet condensate phase in a 
spinor bosonic system |l(j], the stability of spin-ordered 
phases at finite temperature, and unconventional proper¬ 
ties of quantum phase transitions arising from spin cor¬ 
relations P3]. 

As a first step to bridge this gap, here we investi¬ 
gate properties of a strongly correlated spinor Bose gas 
in a 3D optical lattice at both zero and finite tempera¬ 
ture by using a generalized bosonic dynamical mean-field 
theory jl8H2l| . We establish a complete phase diagram 
of lattice spinor bosons with antiferromagnetic interac¬ 
tions, which is shown to contain superfluid (SF), spin- 
singlet condensate (SSC), nematic insulator (NI) and 
spin-singlet insulator (SSI) phases (see Fig.[U). In partic¬ 
ular, we establish the existence of a spin-singlet conden¬ 
sate. At finite temperature, we find a superfluid-nematic- 


insulating phase transition around the tip of Mott lobe 
upon heating, analogous to the Pomeranchuk effect in 
liquid 3 He (see Fig. EH). However, this transition is found 
to be first order for even fillings, but second order for odd 
fillings, contrary to spinless case. 

In sufficiently deep lattices, spin-1 bosons with spin- 
symmetric tunneling are well described by a generalized 
Bose-Hubbard model pj, 

H = —t 'y ' {bl a bj a + bj a bi a ) + —Uq y ' riiirii — 1) 

* 

+ ^/ 2 £(S 2 -2n i )-/r5>,, (1) 

i i 

where b' i(J ( bi a ) is the bosonic creation (annihilation) op¬ 
erator of hyperfine state mp = a at site i , n* = rii a 
with rii C , = b^bta being the number of particles in hy¬ 
perfine state a at site i , Sj = b\ a F aa 'bi a i is the total spin 
operator at site i with F CTCr / being the spin matrices for 
spin-1, g denotes the chemical potential, and t the hop¬ 
ping matrix element between nearest neighbours on the 
lattice. The second term in Eq. m describes a Hubbard 
repulsion between atoms with Uq = (go + 2g 2 )/3, and the 
third term describes on-site spin-dependent interactions 
with U 2 = (g 2 — <7o)/3. Here g s = 4,Trha s /M a with a s 
being the s-wave scattering length in the spin s channel 
and M a the atomic mass. In the following we consider 
the case of antiferromagnetic interaction go < g 2 . This 
is the case, for example, in experiments with 23 Na atoms 

with U 2 /U 0 ^ 0.04 a a. 

In the limit with vanishing spin-dependent interaction 
t/ 2 , the lattice model in Eq. ID reduces to a set of in¬ 
dependent Bose-Hubbard model, where a conventional 
superfluid-insulator transition occurs when Uo/t is tuned 
across the transition point. On the contrary, in the 
limit with t/ 2 Uo, spin-singlet pairs ((S 2 ) = 0) are 
energetically favorable for an even number of atoms in 
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FIG. 1. (Color online) Zero-temperature phase diagram for spin-1 ultracold bosons in a 3D cubic lattice for different anti¬ 
ferromagnetic interactions U 2 /U 0 = 0.01, 0.04 ( 23 Na), 0.3, and 2.0, respectively, obtained via BDMFT (black circle), Gutzwiller 
(red cross) and in Ref. 0 (blue dashed). There are four different phases in these diagrams: superfluid (SF), nematic insulator 
(NI), spin-singlet insulator (SSI) and spin-singlet condensate (SSC). 




each well. In the intermediate, experimentally accessible 
regime with 0 < U 2 — U 0 [8|, [22(, the system has an even 
richer phase diagram as a result of spin-dependent inter¬ 
actions, which compete both with quantum and thermal 
fluctuations resulting in different spin-correlated states. 
We will discuss this regime in more detail in the follow¬ 
ing. 


We investigate the system by means of spinor bosonic 
dynamical mean-field theory (BDMFT), which allows us 
to investigate the system in an unified framework since 
it is non-perturbative and captures the local quantum 
fluctuations exactly. Indeed, BDMFT has been devel¬ 
oped [0| and implemented H-21| successfully for the 
single- and two-component Bose-Hubbard models, where 
it provides a quantitative description for strongly corre¬ 
lated systems in a 3D optical lattice, and the validity of 
this approach has been verified against quantum Monte- 
Carlo simulations [0. Inspired by this success, a spinor 
version of BDMFT is formulated here, where the physics 
on each lattice site is determined from a local effective 
action obtained by integrating out all other degrees of 
freedom in the lattice model Eq. m, excluding the lattice 
site considered. The local effective action is then repre¬ 
sented by an Anderson impurity model, which is solved 
by exact diagonalization [l9(. For comparison, a spinor 
bosonic Gutzwiller mean-field theory is implemented as 
well, and we find remarkable agreement with results from 
the Gutzwiller ansatz in Ref. [15], as shown in Fig. [T] 

Zero temperature — Our main results for the zero- 
temperature case are summarized in Fig. [TJ where a set 


of zero-tenrperature phase diagrams of spinor bosonic 
gases in a cubic optical lattice is shown for different 
anti-ferromagnetic interaction strengths U2/U0 = 0.01, 
0.04 ( 23 Na), 0.3, and 2.0. Four distinct phases, namely 
superfluid (SF), spin-singlet condensate (SSC), nematic 
insulator (NI) and spin-singlet insulator (SSI), are found. 
They are characterized according to the value for the 
condensate order parameter = ( b a ), the nematic or¬ 
der parameter r/> 2 ^ = (S^Ss) — S a p/3(S 2 ), the pair con¬ 
densate order parameter ip\n = (b a bp), and the local 
magnetization M = (S). We remark here that our re¬ 
sults correctly recover unconventional spin ordering both 
in the atomic limit Uq/J = 00 and in the weakly inter¬ 
acting regime. We discuss these phases in detail in the 
following. 

In the deep Mott insulator (MI) with Uq,U 2 t , 
we find that the system favors Mott insulating phases 
with different types of spin order in different interaction 
regimes, i.e. NI characterized by = 0, (j> 2 a p > 0 and 
M = 0; and SSI for an even number of atoms per site as 
a result of the formation of singlet pairs characterized by 
(f>a = 0 and (jy^p = 0 (due to the symmetry of the spin- 
wave function on each site, rii + <S) = (even) 0 , with 
Si being the total spin on site i). Actually, in the strong 
coupling limit Uo,U 2 3> t, one can apply second-order 
perturbation theory and generate an effective spin Hamil¬ 
tonian. Here, for filling N = J2i( n i)/Nut = (odd) where 
Wat is the number of lattice sites, it yields a spin model 
for S = Si/Ni at = 1 bosons 0,0 for describing the 
interplay between ferromagnetic and nematic long-range 
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FIG. 2. (Color online) Zero-temperature SSI-NI-SF (a) and finite-temperature SF-NI-UI phase transitions (b)(c) in a 3D 
optical lattice with chemical potential p/Uo = 1-35 (a)(b) and 0.35 (c), respectively, and interaction U 2 /U 0 = 0.01 (upper left 
panel in Fig. [T] and Fig. Q5J). Inset: zoom of the main figure around the critical point of nematic order (a), and zoom of the 
finite-temperature phase transition (b). BDMFT predicts a first- (second-) order superfluid-nematic-insulating phase transition 
for even (odd) fillings upon increasing temperature, analogous to the Pomeranchuk effect in 3 He. 


order, 

Heff = -Ji £ S 8 • Sj - j 2 £(Si • Sj) 2 , ( 2 ) 

(*i> (ij) 

where (ij) denotes the nearest neighbor sites i, j. One 
can see that the J\ term favors the ferromagnetic phase 
with M 7 ^ 0, but the J 2 term favors the nematic phase 
with (fj^p > 0 and M = 0. The parameter choices 
U 2 /U 0 = 0.04 and U 2 /Uq = 0.3 in Fig. [T] with a total 
filling TV = 1 , 3, for example, correspond to J 2 /J 1 = 1.5, 
1.9 and J 2 /J 1 = 3.3, 3.8, respectively. We find that the 
system is in a nematic insulating phase in both cases, 
which is consistent with the quantum Monte-Carlo pre¬ 
dictions for the above effective spin model 0] . Note that 
the three spin components are not equally populated in 
the nematic phases, while an equal mixture occurs in the 
static mean-field predictions [l]. Similar findings were 
also reported in one dimensional case [oj and could be 
attributed to the quantum fluctuations neglected in the 
static mean-field approach. 

Away from the deep MI regime, where the above re¬ 
sults for Uq,U 2 T§> t based on the perturbation theory are 
questionable, our simulations show that both NI and SSI 
phases are stable against quantum fluctuations within 
the Mott lobes, indicating a higher chance for directly 
observing these phases in realistic experiments at finite 
temperature. For the SF-MI transition, we confirm that 
the transition is first order around the tip of the Mott 
lobe for even filling, while it is second order for odd fill¬ 
ing in a 3D optical lattice jlj, as shown in Fig. Ufa). Note 
that, for small U 2 /Uq (upper-left panel in Fig. [T]), the in¬ 
terplay between tunneling and spin-dependent U 2 yields 
a SSI-NI phase transition before the superfluid-insulator 
transitions occurs 0 . 

Interestingly, we find that upon increasing the spin- 
dependent interaction strength U 2 , the lobe of SSI ex¬ 
pands, while the Mott lobes for odd fillings shrink. Re¬ 


markably, we notice a SSC region, characterized by (fa = 
0, d>lp = 0, (S 2 ) = 0, M = 0 and 0 = (-2 b^h + 6 2 > ^ 
0, which emerges in the phase diagram between the SSI 
Mott lobes after the Mott lobes for odd fillings finally 
vanish when U 2 /Uq > 0.5. Here, we in fact numeri¬ 
cally establish the existence of SSC for spinor bosons in 
a 3D optical lattice (see the right-bottom plot of Fig. [TJ 
where a pronounced region of SSC are found for large 
spin-dependent interaction U 2 /Uq = 2), which was first 
predicted in Ref. Et 

Finally, in the weakly interacting regime with t 
Uq,U 2 , a polar phase with zero magnetization M = 0 is 
found in our simulations, which is consistent with find¬ 
ings in previous works 0 , 0 . 

Finite temperature — One crucial question regarding 
direct observation of the above spin correlated phases 
in realistic experiments is their stability against thermal 
fluctuations. For a typical case with U 2 /Uq = 0.01, the 
influence of thermal fluctuations is shown in Fig. [3] in 
terms of phase diagrams at finite temperatures T/Uq = 
0.005 (left panel) and T/Uq = 0.03 (right panel), where 
five different phases appear in the system, including SF, 
NI, SSI, unordered Mott insulator (UI), and normal state 
(NS) which exists in the low hopping region between the 
insulating lobes and is characterized by fa = 0 but with 
a non-integer filling, as shown by the green filled region in 
Fig. [3] We remark here that these choices of temperature 
are within reach of present cooling schemes, such as spin- 
gradient cooling [25J and the coexistence of these phases 
can be achieved via an external harmonic trap. 

We observe that the SF remains robust against small 
finite temperature, with a tiny change of the insulator- 
superfluid boundary for T/Uq = 0.005. However, for 
higher temperature T/Uq = 0.03, a large shift of the MI- 
SF boundary towards larger hopping is observed. In ad¬ 
dition, upon increasing temperature, we observe a tran¬ 
sition from SF to NI (see Fig. 0(b)), analogous to the 
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FIG. 3. (Color online) Finite-temperature phase diagram for spin-1 ultracold bosons in a 3D cubic lattice with an anti¬ 
ferromagnetic interaction U 2 /U 0 = 0.01 for different temperatures T/Uo = 0.005 (left) and T/Uo = 0.03 (right). There are five 
different phases in this diagram: superfluid (SF), nematic insulator (NI), spin-singlet insulator (SSI), unordered insulator (UI) 
and normal state (NS). SSI is fragile against temperature and shrinks into NI with increasing temperature, while UI develops 
from the lower hopping regions of NI, due to thermal fluctuations. 


Pomeranchuk effect in liquid 3 He where the system fa¬ 
vors localization upon heating at low temperature, since 
spin fluctuations in the Mott insulating phase can carry 
more entropy than the superfluid. Upon further increas¬ 
ing temperature, a NI to UI transition is observed. In¬ 
terestingly, we find that the superfluid-Mott-insulating 
phase transition for even fillings is first order, as shown 
in Fig. (3](b). The physical reason is that the Mott state 
with even filling prefers the lowest spin state (spin 5 = 0), 
while the SF contains higher spin states, similar to the 
zero-temperature case [Fig. dja)]. 


For the SSI phase (N = 2 with U 2 /U 0 = 0.01), we find 
that it is sensitive to finite temperature, and the under¬ 
lying physics is that the spin-gap U 2 term is small and 
easily destroyed by thermal fluctuations. As shown in 
the right panel of Fig. [31 for example, no SSI occurs at 
temperature T/Uo = 0.03. For experiments with 23 Na 
(U 2 /U 0 ~ 0.04) in a 3D cubic lattice formed by laser 
beams of wavelength 1064 nm and intensity Vj ss 16 Ur, 
we find the critical temperature to be T c s=s 1 nK, where 
Ur denotes the recoil energy. In addition, the nematic 
phase, NI (the critical temperature is T c ss 5 nK, for the 
same experimental setup with 23 Na but with an inten¬ 
sity of Vi ~ 14Ur), is reduced with increasing temper¬ 
ature, in favor of developing a non-ordered Mott state. 
We remark here that these correlated insulating states 
can be detected from their excitation spectra or density 
correlations 12, 1(|, by using, for instance, Bragg scat¬ 
tering Q,. quantum gas microscopy [26[ or optical bire¬ 
fringence [16|, as long as the lifetimes of these states of 
spinor lattice bosons are sufficiently long Q. In recent 
experiments, a long lifetime steady state of spinor bosons 
in optical lattices was reported §, as well as the existence 
of spin-nematic ordering in a spherical trap [27 ]. 


In conclusion, we have investigated quantum phases of 


ultracold spinor Bose gases loaded into a cubic optical 
lattice, based on spinor bosonic DMFT. We obtain the 
complete phase diagram of spinor bosons with antiferro¬ 
magnetic interactions in an optical lattice at both zero 
and finite temperature. In particular, we demonstrate 
the existence of a spin-singlet condensate. Interestingly, 
we observe that the superfluid can be heated into a Mott 
insulator with even (odd) filling via a first- (second-) or¬ 
der phase transition, analogous to the Pomeranchuk ef¬ 
fect in 3 He. We find that the critical temperature of 
ordered states, such as nematic phase and spin-singlet 
insulator, is within reach of present cooling schemes, in¬ 
dicating the chance to directly observe these phases using 
current experimental techniques. 
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